THE WILSON FUNCTION TRANSFORM 



WOLTER GROENEVELT 

Abstract. Two unitary integral transforms with a very- well poised 7_F6-function as a kernel are 
given. For both integral transforms the inverse is the same as the original transform after an 
involution on the parameters. The 7_F6-function involved can be considered as a non-polynomial 
extension of the Wilson polynomial, and is therefore called a Wilson function. The two integral 
transforms are called a Wilson function transform of type I and type II. Furthermore, a few explicit 
transformations of hypergeometric functions are calculated, and it is shown that the Wilson function 
transform of type I maps a basis of orthogonal polynomials onto a similar basis of polynomials. 



1. Introduction 

The Jacobi polynomials are orthogonal polynomials that that are solutions to the hypergeo- 
metric second order differential equation. The Jacobi polynomials have an explicit expression as 
2-Fi-hypergeometric series, see e.g. [Q. Wilson |23| gave a generalization of the Jacobi polynomials 
as orthogonal 4i ? 3-polynomials. These polynomials are nowadays called the Wilson polynomials. 
The Wilson polynomials no longer satisfy a second order differential equation, but a second order 
difference equation. Another way to generalize the Jacobi polynomials is to consider non-polynomial 
solutions of the hypergeometric differential equation. Spectral analysis of the hypergeometric dif- 
ferential equation leads to a unitary integral transform, called the Jacobi function transform, see 
e.g. |12j . The kernel in this integral transform is given by a non-polynomial 2 -^-functions called 
the Jacobi function. Koornwinder mentions in [111 §9] that there probably also exist functions, 
naturally called Wilson functions, that generalize the Wilson polynomials as well as the Jacobi 
functions. 

In this paper we consider non-polynomial eigenfunctions to the Wilson second order difference 
operator. To study the associated Wilson polynomials, Ismail et al. 5 , and Masson ^1], showed 
that general solutions of the difference equation are given by very-well-poised 7-Fg-functions. The 
analytic part of a certain non-polynomial solution we call a Wilson function. We show that the 
Wilson functions are the kernel in two unitary integral transforms, which we call the Wilson function 
transform of type I and type II. The Wilson function satisfies the property, called the duality 
property, that after an involution on the parameters, the geometric and the spectral parameter 
can be interchanged. The parameters obtained from the involution are called the dual parameters. 
The inverse of the Wilson function transform of type I, respectively type II, is again a Wilson 
function transform of type I, respectively type II, with dual parameters. The transforms can be 
made completely self-dual by choosing a fixed point of the involution. 

The Wilson function transform of type I is a generalized Fourier transform with respect to the 
same measure as the Wilson polynomials. This measure consists of an absolutely continuous part 
supported on [0, oo), and a finite number of discrete mass points. The Wilson function transform of 
type II is a generalized Fourier transform with respect to a one-parameter family of measures which 
consist of an absolutely continuous part supported on [0,oo), and an infinite number of discrete 
mass points. The extra parameter labels the different sets of discrete mass points. 

The Wilson functions can be considered as a formal limit case of the Askey- Wilson functions for 
q — > 1. The Askey- Wilson functions are eigenfunction to the Askey- Wilson second order difference 
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operator for < q < 1. The Askey Wilson functions are the kernel in the Askey- Wilson function 
transform, which is found by Koelink and Stokman j^j. They show in JO] that the Askey- Wilson 
function has an interpretation as spherical function for the quantum SU(1, 1) group. From this 
point of the view the Askey- Wilson function is a g-analogue of the Jacobi function, since the Jacobi 
function has similar interpretation for the Lie group SU (1,1), see |12j . So the Wilson functions 
in this paper give a new limit case of the Askey- Wilson functions. We do not know if the Wilson 
function also has an interpretation as a spherical function. 

The Wilson function can also be considered as a formal limit case of Ruijsenaars' i?-function 
[T7j . [T8] . [T9] . which is an eigenfunction to the Askey- Wilson second order difference operator for 
\q\ = 1. The i?-function is given by a Barnes- type integral which is considered as a generalization of 
the Barnes integral representation for the 2-^i-series. Using the Barnes-type integral representation 
for the 7-F6-series |2U1 (4.7.1.3)], the i?-function can also be considered as a generalization of the 
Wilson function. 

In a future paper we will show that the Wilson functions, and also the Wilson polynomials, 
have an interpretation as Racah coefficients for tensor products of positive discrete series, negative 
discrete series, and principal unitary series representations of the Lie algebra su(l, 1). Both Wilson 
function transforms in this paper have an interpretation in the context of Racah coefficients. 

The organization of this paper is as follows. In section [21 we give some well-known properties of 
the Wilson polynomials. The Wilson polynomials are eigenfunctions of a second order difference 
operator A, and we show that the Wilson polynomials are also eigenfunction of the same difference 
operator with dual parameters. 

In section |21 we consider a certain type of non-polynomial eigenfunction of a difference operator 
L which is closely related to A. These eigenfunctions, the Wilson functions, are also eigenfunctions 
of the difference operator L with dual parameters. 

In section HI we define a Hilbert space A4, and, using the asymptotic behaviour of the Wilson 
function, we show that a truncated inner product of two Wilson functions approximates a reproduc- 
ing kernel. This leads to a unitary integral transform, which we call the Wilson function transform 
of type I. 

In section we define a different Hilbert space 7i. Again using asymptotic behaviour of the 
Wilson function, we show that a truncated inner product of two Wilson functions approximates a 
reproducing kernel. This leads to the Wilson function transform of type II. 

In section El the Wilson function transforms of a Jacobi function, and of a Wilson polynomial, 
are calculated explicitly, using an integral representation of the Wilson function. Also we show that 
the Wilson function transform of type I maps an orthogonal system of Wilson polynomials onto 
the same orthogonal system with dual parameters. 

Notations. We use the standard notation for the hyper geometric series, i.e. 



A hypergeometric series is called very- well poised if p = q + 1, ai + 1 = b\ + 02 = . • • = b g + a q+ \ 
and 61 = 02/2. For a very-well poised 7-Fg-series of argument 1 we use Bailey's W notation, see 
i.e. 




where (a) n denotes the Pochhammer symbol, defined by 



. , T(a + n) . . . , 
(a)„ = y ' = a(a + l)(a + 2) . . . (a + n - 1) 




n € Z>q. 



W(a; b, c, d, e, 



f) = rF 6 



a, 1 + ia, b, c, d, e, f 
|a, 1 + a — b, 1 + a — c, 1 + a — d, 1 + a — e, 1 + a — f 



;1 
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If the series does not terminate, the condition for convergence is 5R(2a + 2 — b — c — d — e — f) > 0. 
Also, since the product T(a + b)T(a — b) frequently occurs in this paper, we use for this product 
the shorthand notation Y(a ± b). 

2. Wilson polynomials 

In this section we recall some well-known properties of the Wilson polynomials. The Wilson 
polynomial satisfies a second order difference equation in its degree (the three-term recurrence re- 
lation), and also a second order difference equation in its argument. The goal of this section is to 
point out that the two difference equations are the same after a change of the parameters. 

The Wilson polynomials R n (x), see |23j . §6.10], [HJ, are polynomials in x 2 of degree n. They 
can be defined by the initial values R-i(x) = 0, Ro(x) = 1, and the recurrence relation 

-(a 2 + x 2 )R n (x) = C n [R n+1 (x) - R n (x)] + D n [R n ^(x) - R n (x)] , (2.1) 

where 

[n + a + b + c + d — l)(n + a + b) (n + a + c) (n + a + d) 



C n = C n (a,b,c,d) 
D n = D n (a,b,c,d) 



(2n + a + b + c + d - l)(2n + a + b + c + d) 
n(n + b + c - l)(n + b + d - l)(n + c + d - 1) 



(2n + a + b + c + d - 2)(2n + a + 6 + c + d - 1) 

The explicit expression for the polynomials i? n is given by 

/-n.n + alHcld-M + K.a-K \ 
-RnW = Rn{x;a,b,c,d) = 4 F 3 ;1 . (2.2) 

\ a + b,a + c,a + d J 

Let a, 6, c, d € C be such that non-real parameters appear in conjugate pairs with positive real 
part, and such that the pairwise sum of any two parameters has positive real part. We define the 
measure dp,(x) = dfi(x; a, b, c, d) by 



/I roo _^ 
f(x)dfi(x) = — f(x)w{x)dx + i^2f{ 

w(x) = w(x; a, b, c, d) 



j ™j T 

where 

r(a ± ix)T(b ± ix)T(c ± ix)T{d ± ix) 
T(±2ix) ' 

The points Xk are of the form xj~ = i(e + k), where e is any of the parameters a, b, c, d with e < 0. 
The sum is over k € Z>o, such that e + < 0. The weights Wk are the residue at x = Xk of w(x). 
In particular, if all parameters are positive, or occur in pairs of complex conjugates with positive 
real part, the measure d[i is absolutely continuous. The Wilson polynomials are orthogonal with 
respect to the measure dfx, i.e. 

R m (x)R n (x)dfJ,(x) = 

a + b + c + d-1 n\T(a + b)T(a + c)T(a + d)T(b + c + n)F(b + d + n)T(c + d + n) 



5 m 



2n + a + b + c + d-l (a + b) n (a + c) n {a + d) n (a + b + c + d - l) n T{a + b + c + d) 

(2.3) 



The Wilson polynomials also satisfy a difference equation in x, given by 

n(n + a + b + c + d- l)R n (x) = A(-x) [R n (x + i) - R n (x)] + A(x) [R n (x - i) - R n (x)] , (2.4) 
where 

, , (a + ix)(b + ix)(c + ix)(d + ix) 
^ ~ 2ix(2ix + 1) ' 
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We define the difference operator A by 

A = A(x)(T„ i -I)+A(-x)(T i -I), (2.5) 

where I denotes the identity operator and T is the shift operator (i.e. T z f(x) = f(x + z)). From 
the difference equation (|2.4|) it follows that the polynomials R n (x; a, b, c, d) are eigenfunctions of A 
for eigenvalue n(n + a + b + c + d — 1). 

The recurrence relation can be written in a self-dual way. Given the parameters a, b,c,d € C, 
we define dual parameters by 

1 1 

a = -(a + b + c + d- 1), b = -(a + b - c - d + 1), 

I . i ( 2 - 6 ) 

c = -(a - b + c - d + 1), d = -(a-b-c + d + l). 

It is an easy verification that (a, 6, c, d) i— > (a, 6, c, d) defines an involution on the parameters. For 
a function / = f(a,b,c,d), we define / by / = f(a,b,c,d). We use the same notation for other 
objects, like measures, sets and operators. Denote the Wilson polynomial R n {x) by P\(x), where 
A = i(n + a), i.e. 

_ , . /a + iA, a — iA, a + ix, a — ix N 

-Pa (a;) =4^3 , , , ; 1 



a + b,a + c,a + d 

then we see that P\(x) = Pc(A), since a + e = a + e, for e = b, c, d. The recurrence relation (|2,lj) 
and the difference equation (|2.4j) for the Wilson polynomials can now be written as 

AP x (X) = -(a 2 + x 2 )P x (\), 

AP x (x) = -(~a 2 + X 2 )P x (x). 

So we see that, for A G i(a + Z>o), the Wilson polynomial P\(x) is an eigenfunction of both A and 



A. 



3. Wilson functions 



In the previous section we observed that, for A 2 = — (d+n) 2 , the Wilson polynomials are solutions 
to the eigenvalue equation 

(Af)(x) = -(\ 2 + a 2 )f(x), (3.1) 

For more general values of A, solutions to (|3.1|) can be given in terms of very-well poised -jFq- 
series. This is shown by Ismail et al. and by Masson who investigate the associated Wilson 
polynomials. Let ip\{x) be the function defined by 

r(6 + c)r(5 + b + c + iA)r(l - d + »A)r(l -d±ix) 



ijj x (x) = il>\(x;a,b,c,d) 



T(b + c + i\±ix) 
W(a + b + c — 1 + i\; a + ix,a — ix,a + i\,b + i\,c + iA). 



The dual parameters a,b,c,d are still defined by 1)2.6)) . By [SJ Thm. 2], or (2.5)], the function 
4>\(x) is a solution to the eigenvalue equation 1)3.1)1 . 
Instead of t^a (^) > we study the closely related function 

a r \ a i u a\ T(a + ~b + c + i\) 

(px{x) = 4>x(x;a,b,c,d) - 



T(a + b)T(a + c)T(l + a - d)T(l -d- i\)T(b + c + i\±ix) (3.2) 
x W{d + b + c — 1 + i\;a + ix, a — ix,d + i\,b + iA, c + iX). 
So <j>\(x) = tp\(x)/ K(x, A), where X(x, A) is the function given by 

K(x, A) = T(a + 6)r(a + c)r(6 + c)T(l + a - d)r(l - d ± ix)T(l - d ± iA). 
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We call the function <p\{x) a Wilson function. The Wilson function 4>\(x) has the advantage that 
it is symmetric in a, b, c, 1 — d, cf. B,emark I4.5f iii). while ip\(x) is not. Also 4>\(x) is an analytic 
function in (x, A) G C 2 . 

The 7i*6-series in the definition of (f>\(x) converges absolutely for 5J(1 — d — iX) > 0. Writing 
the 7-Fg-series as a sum of two balanced 4F3-series, we have an expression for 4>\(x) which always 
converges: 

r(l — a — d) fa + ix, a — ix, a + iX, a — iX 

(p\{x) = = 4-F3 , ;i 

T{a + b)T{a + c)T{l-d±ix)T{l-d±iX) \ a + b,a + c,a + d 



T(a + d-l) 



(3.3) 



x 4^3 



r(l + b - d)T(l + c - d)T(a ± ix)T(d ± iX) 

l-d + ix,l-d-ix,l-d + iX,l-d-iX _\ 
1 + b-d, l + c-d,2-a-d ' J ' 

This follows from 2, §4.4(4)] with parameters specified by 

a 1— > a + 6 + c — 1 + iX, c^c + iX, d^b + iX, e 1— >a + iA, f^a — ix, g\-^a + ix. 

From H3.3|) we see that <^a(^) is an analytic function in (x, A) G C 2 . Observe that for A = ±i(a + n), 
the second term in (|3.3j) vanishes because of the factor T(d ± iA) -1 , and then we see that V , a(^) = 
iC(x, A)^>a(^) reduces to a Wilson polynomial. So tp\(x) is the analytic continuation of the Wilson 
polynomial in its degree. 

From (|3.3|) and (|2.6|) follows the duality property 

4x{x) = 4> X (X). (3.4) 

This duality property is similar to the duality property for the Askey- Wilson functions in For 
the 7-F6-series in the definition of the Wilson function, ()3.4|) is implied by Bailey's transformation 
§7.5(1)]. 

Since <fi\(x) = ip\(x)/ K(x, A) and ip\(x) is a solution to eigenvalue equation (|3.1|) . the Wilson 
function satisfies the equation 

M kIx,x) A M K(x,x) <Px{x) = -(A 2 + a 2 )0 A (x). 

Here M K ^ X ^) denotes multiplication by K(x,X). From this we obtain the following proposition. 

Proposition 3.1. The Wilson function (fi\(x) is a solution to 

(L/)(x) = (A 2 + a 2 )/(x), 

where L is the difference operator defined by 

L = B{-x)Ti + [A{-x) + A(x)] I + B{x)T„i, 

(a + ix)(b + ix)(c + ix)(d + ix) 



A(x) 
B(x) 



2ix(2ix + 1) 
(a + ix)(b + ix)(c + ix){l — d + ix) 



2ix(2ix + 1) 

In the next two sections we consider the action of the second order difference operator L on two 
different Hilbert spaces, which leads to two different integral transforms with the Wilson function 
as a kernel. The method we use comes down to approximating with the Fourier transform, using 
asymptotic expansions of the Wilson functions. This method is essentially the same method as 
used by Gotze [Bj, and Braaksma and Meulenbeld [3], for the Jacobi function transform. A similar 
method is also used by Koelink and Stokman in j^j for the Askey- Wilson function transform, and 
in [7] for the continuous Hahn function transform. 
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4. The Wilson function transform: type I 

4.1. The Hilbert space .M. Let V C C 4 be the set of parameters (a, 6, c, d) satisfying the following 
two conditions: 

• The parameters a, b, c, 1 — d are real except for pairs of complex conjugates with positive 
real part. 

• The pairwise sum of a, b, c, 1 — d is contained in C \ (— oo, 0]. 

A direct verification shows that the assignment (a,b,c,d) h- > (a,b,c,d), see (J2.6[) . is an involution 
on V. Throughout this section we assume (a, 6, c, d) G V. 
Let M be the weight given by 

, ,, . , ,. , , N T(a ± ix)T(b ± ix)T(c ± ix)T(l — d ± ix) 
M(x) = M(x;a, b, c, d) = -± ^ T(±2ix) ' 

The weight M is positive for Observe that M(«; a, b, c, d) = w(-; a, b,c,l — d), where w is the 

weight function for the Wilson polynomials. We assume that the function M has only simple poles. 
This imposes conditions on the parameters a, b, c, d that can be removed afterwards by continuity 
in the parameters. For e G C define the set T> e by 

V e = {i(e + n) | n G Z>o, e + n < 0}, 

and let P = V a U Df, U D c U T>i-d- We define the measure dm(-) = dm(-;a, b, c, d) by 

/ /(x)dm(x) = - f /(x)M(x)dx + iV /(x)ResM(z). 
J 2lT JO z=x 

If x € V a , we have explicitly 



xev 



. . r(a + 6)r(a + c)r(l + a - d)r(6 - o)r(c - a)r(l - d - a) 
i Res M{z) = — 

^=j(a+n) 1 ( — zaj 

(2a) n (a + l)n(a + &)n(a + c)„(l + a - d) n 



n\ (a) n (l + a - 6) n (l + a - c) n (a + d) n 

and then we see that for (a, 6, c, d) € V, the measure dm is positive. Recall that (f>\(x) is symmetric 
in a,b,c, 1 — d, and by (|3.3|) the Wilson function is obviously even in x and A, therefore (j>\{x) = 
<frv(x)- From this it follows that for x G supp dm(-) and A G supp dm(-) the Wilson function 4>\(x) 
is real valued. 

We define the Hilbert space M. = A4(a, b, c, d) to be the Hilbert space consisting of even functions 
that have finite norm with respect to the inner product (•, -)m defined by 



(/, 9)M = j f(x)g(x)dm(x). 
4.2. The Wronskian. For < N < oo, we define a pairing (•, by 

1 r N 

(f,g)N = 7T / f(x)g(x)M(x)dx + i V f(x)g(x)Re S M(z). 

If / and g are real valued functions in M, the limit TV — > oo gives the inner product (f,g)M- 
For functions /, g that are analytic in C, we define the Wronskian [f,g] by 

[/, <7]0) = ^ - i) - /(x - i)g(x)} B(x)M(x)dx. 

Lemma 4.1. Let f,g be analytic in C and even, then z i— > [f,g](z) is odd in z. 
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Proof. Let I be the function given by 

I(x) = (f(x)g(x - i) - f(x - i)g{x))B{x)M{x), 

then [f,g](z) = f* +l I(x)dx. Since f(x),g(x) and M(x) are even functions in x, and B(—x)M{x) = 
B(x + i)M(x + i), we have I(—x) = —I(x + i). Therefore 

rz+i r—z~i p — z—i p—z+i 

I I(x)dx = — I I(—x)dx= / I(x + i)dx = — / I(x)dx. 

Jz J—z J—z J—z 

Hence z \— > [/, </](#) is an odd function in z. □ 

Proposition 4.2. For N and /or euen analytic functions f and g we have 

{Lf,g) N -(f,Lg) N = [f,g](N). 

Proof. Recall that we assume that the poles of M(x) are simple. For even functions / and g we 
have 

(f,g)N = ^~ [ f(x)g(x)M(x)dx, 
47r Jc N 

where Cat is a contour in the complex plane defined as follows: 

• Cn starts at x = —N and ends at x = N, 

• Cn is invariant under reflection in the origin, 

• Cn separates the sequence of poles i(a + n), n G Z>o from the sequence —i(a+m), m S Z>o, 
and similarly for poles of M(x) corresponding to b, c, 1 — d. 

Now we have 

' 1 (Lf)(x)g(x)M(x)dx - f f{x){Lg)(x)M(x)dx = 



4vr Je N ^ Jc N 

j- [ (f(x + i)g(?) - f(x)g(x + »)) B{-x)M(x)dx 
J Cm 

j c (f( x ~ i)9{x) ~ f( x )g(x ~ «)) B(x)M(x)dx. 
Using B{—x + i)M(x — i) = B(x)M(x), we can write this as 

^( I - I ) (f( x ~ i)9{x) - f(x)g(x - tj) B(x)M(x)dx 



Cn JCfj+i 



Here Cn + i = {x £ C\x — i £ Cn}- The integrand has its poles at x = i(e + 1 + re), x = —i(e + m), 
for re, m G Z>o and e = a, b, c, 1 — d. Now we make a closed contour by connecting Cn and Cn + « 
at the end points in a straight line (there are no poles of M on these lines for iV large enough), 
then the integrand I(x) has no poles inside the closed contour. So, by Cauchy's Theorem, 

-N+i r N+i 



and from this we obtain 



Cn JC N +i J-N JN 



(Lf,g) N - (f,Lg) N = \[f,g](N) - \[f,g](-N). 

Now the proposition follows from Lemma 14.11 □ 

Since the Wilson functions are eigenfunctions of L for eigenvalue a 2 +A 2 , we find from Proposition 
.21 the following. 
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Proposition 4.3. For X ^ X' , 

[Mxr](N) 
((pxAxON = X 2_ X ,2 ■ 

Next we want to let N — > oo in Proposition 14.31 so we need the asymptotic behaviour of the 
Wilson function and of B(x + iy)M{x + iy) for x — > oo and < y < 1. We have 

M(x + iy) = l 67r 3 a .2a + 2 6+ 2c-2d-l e -2 7 rx-2i 7 r 2/ A + o(-)) , 
2/ ,.x V 



B(x + iy) = -—[l + 0[- 
4 \ Vx 

The asymptotic behaviour of the weight function M can be obtained from |16l §4.5] 

£|±|l = ^-^l + -L(a-6)(a + 6-l) + 0(l)), |*| -+ oo, |arg(z)|<vr, (4.2) 

and from applying Euler's reflection formula for the T-function 

7rr(a + ix) 



r(a ± ix) 



r(l — a + ix) sin7r(a — ix) 



To determine the asymptotic behaviour of 4>\(x), we expand 4>\(x) in functions with nice asymptotic 
behaviour. 

Proposition 4.4. 

<£ A (x) = c(A)$ A (x) + c(-A)$_ A (x), 

where 

_ 1 / 5 + iA, 6 + iA, c + iA, 1 - J + iA \ 

A ~ r(6 + c + iA±ix) 4 3 \^ + c + iA + ix, b + c + iA - ix, 1 + 2iA ' J ' 

c(A) - r <- 2!A) 



r(2 - iA)r(6 - i\)T(c - iX)T(l -d-iX) 

Proof. This follows from transforming the yi^-function in the definition Q3.2JI of (f) x (X) by [2J §4.4(4)] 
with parameters specified by 

a i— > a + 6 + c — 1 + iA, c i— > a + ix, d i— > a — ix, e^a + iX, ft— >6 + iA, gi-^c + iA. 

□ 

Remark 4.5. (i) Note that the dual weight function M can be expressed in terms of the c-function 

as 

(ii) The functions <1?a an d ^-a are i n general not solutions to the eigenvalue equation in Propo- 
sitioning Let us define 
* A (x) = 

r(l - a - ix)T(l -b- ix)r(l - c - ix)T(l -d- ix)T{2 - d - ix + 2iA) 

T(l + a - d + iX - ix)T(l + d- d + iX- ix)T(2 -b-d + iX- ix)T(2 -c-d + iX- ix)T(d - ix) 

sin7r(l + a — d + ix + iX) TTr/ _ , „ ,.„_.,., 7 ., r _ ., N 

x — r W (1 — a — zx + 2zA; 1 — a — zx, 1 — a + zA, 1 — a + zA, + iX, c + iX). 

sm7r(a — ix) 
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From |141 (2.12)] it follows that and *&-\ are solutions to the eigenvalue equation (jH.ljl . so 
ty\(x)/K(x, A) is an eigenfunction of L for eigenvalue A 2 + a 2 . The Wilson function (j>\{x) can be 
expanded in terms of ^>\{x) and *$>-\(x) as follows: 

r(l - d± ix)4>x{x) = c(X)V x (x) + c(-A)*_ A (x), 

where the c-function is the same as in Proposition 14.41 This expansion follows from |2(Jl (4.3.7.8)], 
with parameters given by 

a i — ► 1 — d ~\~ 2iX — ix, b i— > 1 — d — ix, c *— > 1 — a + iX, 
d^l — d + iX, e i— > 6 + i A, / i — > c + zA, 

and Euler's reflection formula for the T-function. From this expansion is it possible to determine the 
asymptotic behaviour of <j)\(x) for x — > oo. We will use the simpler functions <£±a from Proposition 
14.41 to compute the asymptotic behaviour of the Wilson function. 

(hi) Writing c = \{d — b + c — d + 1) in Proposition 14.41 we see that the Wilson function 
4>\(x; a, b, c, d) is symmetric in a, b, 6,1 — d, and by the duality property (|3.4|) then also in a, b, c,l—d. 

To determine the asymptotic behaviour of &\(x) we use Euler's reflection formula to rewrite the 
T-functions in front of the 4i ? 3-function and we use (|4.2|) . Then we find for <E>a, f° r V £ K an d 

x — ► oo, 



$ A ( X + iy) = J_ a .d-a-6-^MA eTO -M W1 / ( x + J/_ ( a + 6 + c _ d + 2U ) + Q( J_ ) j . (4.3) 
z^7^ \ 2X x z ' ' 

From Proposition 14.41 it follows that the asymptotic behaviour of the Wronskian [<ft\,<ft\'](N) can 
be computed from the four Wronskians [&±\, $±\'](N). So we must compute [$>\, $>\r](N) for 
N -> 00. 



Lemma 4.6. For N — > 00, 



[$ A , * A ,](JV) = »(A - A')iV- 2i(A+A ^ ( 



l + O 



Proof. Let (?(x) be the function given by 

G(x) = <S>x{x)$\'{x ~ i) ~ - %)®\>{x). 
From the asymptotic behaviour of <3?a given above we find, for < y < 1 and x — > 00, 

G(x + iy) = - (A 2 ~ 2 A/) x 2d-2a-2b-2c-l-2i(X+X') e 2^+2iny A + _ 

Using the asymptotic behaviour of i?(x + iy)M(x + %) for x — > 00 gives 

G(x + iy)B{x + iy)M(x + iy) = 2tt(X - A')2r 2i(A+A) M + o{^j V 

Note that the main term in de asymptotic expansion is independent of y. Next we write 

* Al $ A ' PO = 77- / G(x)£?(x)M(x)dx = -7T-- G(N + iy)B(N + iy)M(N + iy)dy, 
2vr J N 2m J 

and we apply dominated convergence to obtain the result. □ 
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4.3. Continuous spectrum. In this subsection we assume A, A' 6 R, and since 4>\ is even in A, 
we may assume A, A' > 0. 

Proposition 4.7. Let f be a continuous function, satisfying 



/(A) 
Then 



0(X 5 ), AjO, <5>0. 



lim -L f°° f(\)(4>\, 4>\')n dX — • /(Y) 



N^oo2irJ Q M(X') 
Proof. Prom the c-function expansion in Proposition 14.41 we find 

[4 >x , ( f )xl }(N)= Yl c(eX)c^X')[^ eX ,^y](N), 

e^6{-l,l} 

and then we obtain from Lemma 14.61 and Proposition 14.31 



\<P\,<P\>)N = « 

e,fe{-i,i} 

We multiply both sides with an arbitrary function /(A), and we integrate over A from to oo. The 
function / must satisfy certain conditions that we determine later on. Letting N — > oo then gives 

;>oo 

lim / f(X)((j>\,(t>\')Nd\ = 

N^OO In 



lim i 



where 



/(A){^i(A) cos (2[A + A'] ln(AO) + ^j(A) sin (2[A + A'] ln(iV)) 

+ *(A) cos (2[A - A'] In(iV)) + ^(A) ^ ^\~_ X '^ } dX 

^ = aTa 7 ( a(A)c ~ (y) " c ~(-a)c(-a')), 

"^2 (A) = ^ (c(A)c(A') + c(-A)c(-A0) , 
^3 (A) = — ^ (c(A)c(-A') - c(-A)e(A')), 

^ 4 (A) = -i(c(A)c(-A') + c(-A)c(A')). 

Observe that c(A')c(— A') — c(— A')c(A') = 0, so ^3 has a removable singularity at A = A'. From 
the Riemann-Lebesgue lemma we find that the terms with ipi, i = 1,2,3, vanish, provided that 
fipi G L 1 (0, 00). We recognize the term with ip^ as a Dirichlet integral. Using the well-known 
property (see e.g. |221 §9.7]) for Dirichlet integrals 

lim - f° g(x f n[t{x - V)] dx = g(y), (4.4) 
t^oo 7T J x-y 

for a continuous function g G L 1 (0, 00), we obtain 



DC 



lim / f(X)(<t> x ,ct>x} N dX = ^4(A')/(A') = 2vrc(A')c(-A')/(A / ) = 2nS^- 

N->oo_l Q M(X') 
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The asymptotic behaviour of tpi, for i = 1,2, 3, 4, can be obtained in the same way as the asymptotic 
behaviour of M, see (|4.1j) . Then we find, for i = 1, 2, 3, 4, 



l^(A)| 



0(A-2- s -^+V A ), A -» oo, 
©(A- 1 ), A |0. 



So, if / satisfies the conditions given in the proposition, then ftjji E -L 1 (0, oo). □ 
Let / be a continuous function. We define a linear operator T by 

(^/)(A) = / /(x)^(x)dm(is). 



We call T the Wilson function transform of type I. We denote the continuous part of the above 
integral by T c f , i.e. 



oo 



o 



From the asymptotic behaviour of and M(x) we find that both J-f and .T^/ are well defined 
if / satisfies the conditions 



/(*) 



p^rf-a-t-c-egB^ x^OO, e > 0, 

Oix 5 - 1 ), x I 0, 5 > 0. 



Observe that if / satisfies the condition given above, then / 6 j\4. Let .Mo C M be the space 
consisting of continuous functions satisfying the above conditions. Then Mq is a dense subspace of 
A4 (it contains for instance the Wilson polynomials R n (x; a, b,c,l — d), which form an orthogonal 
basis for M). 

Proposition 4.8. Let g G M and A > 0, then (^(i^?)) (A) = g(X). 

Proof. Define /(A) = M(X)g(X), then / satisfies the conditions given in Proposition 14.71 Then we 
have 

[T{F c g))(\') = j Mx){^ f°° g{\)^{x)M{\)d\) dm(x) 



o 



lim — / /(A) — / <f> x (x)(l)x(x)M(x)dx + i V A (x)0 v (x)ResM(z) dA 



i r 00 

lim — / f(\)((j)\,<j>x>)Nd\ 



N 

/(A' 



M(A' 



Note that the first integral converges absolutely for <? € A4q, so in that case interchanging the order 
of integration is allowed. □ 



In the next subsection we show that the dual Wilson function transform T is the inverse of the 
Wilson function transform T. To do this, we must consider the discrete spectrum of the difference 
operator L. 
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4.4. Discrete spectrum. From the asymptotic behaviour of &\(x) and M{x), see (|4.3jl and (|4.1|l . 
we obtain 

|$a(^)| 2 M(x) =C(x 43(A) " 1 ), rr^cxo. 

So for 3(A) < we have <&\ E A4. In this subsection we assume that A E V and that the set V is 
not empty, so 9(A) < 0. In this case c(— A) = 0, and therefore 4>\{x) = c(A)<3? A (x) E A4. First we 
show that 0a is orthogonal to </> A ' if A' 7^ A. 

Proposition 4.9. For A E T> , A' E supp(dm), A ^ A' ; we /iaue 

(</>A, 0A')A4 = 0- 
Proof. From Propositions 14.31 and 14.41 we obtain 

c(A)c(AQ[cD A ,cD A ,](iV) + c(A)c(-AO[cD A ,cD_ A ,](iV) 

{4>\,4>\')n = A 2 _ A / 2 • 

Then Lemma 14.61 gives for large N 

i{\ - \')c{X)c{\')N- 2i{ ~ x+x '^ + i(X + A')c(A)c(-A , )A f - 2i{A - A ' ) 



(A + A0(A-A') 



l + O 



(4.5) 



Recall that for A E £> we have A E iM<o. Then it is clear that for A'eR the right hand side tends 
to zero for N — > 00. In case A' E D the second term vanishes, and we have ^s(X + A') < 0. So in 
this case the right hand side also tends to zero for — > 00. □ 

It remains to calculate the squared norm of (f>\ in case A E T>. 

Proposition 4.10. For A E V 

{^xA\)m = (i Res M(A') 



A'=A 

Proof. We use expression (|4.5j) . where we let A' — > A. Then for large N 

I AA\ /„ 1 \ _1 



lim 

A'^A 



Letting N 



A, 4>\')n : 

00 gives the result. 



2A c(A,c(A)^ +i£ (A,(R S . ( _ A , ) 



l + O 



N 



□ 



Remark 4.11. If A, A' E 2?a then Propositions 14.91 and l4~10l give orthogonality relations for a fi- 
nite number of Wilson polynomials with respect a measure that has only finitely many moments. 
Explicitly, if we assume a, b, c, 1 — d > 0, we have for n, m < ^(1 — a — b — c — d) 

2 

dx = 



— / R „(•'• )./-?,»(■'•) 



T(a + zj;)r(6 + ix)T(c + ix) 



T(l - d + ix)r(2ix) 

a + b + c + d- 1 n! (b + c) w (fe + d) n (c + rf) w 

nm a + 6 + c + d+2n-l(a + b) n (a + c) n (a + cf)„(a + b + c + d - l) n 

T(a + 6)r(a + c)T(b + c)r(l -a-6-c-d) 
X r(l - a - d)T(l -b- d)r(l - c - d) ' 

where R n is the Wilson polynomial defined by (|2.2j) . Neretin [E] found this orthogonality relation 
using an explicit evaluation of the Barnes-type integral 



1 

2^ 



r(a + ix)T(b + ix)T{c + ix) 



dx. 



T(l-d + ix)T(2ix) 

Note that the above orthogonality relations remain valid for n + m < 1 — a — b — c 



d. 
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4.5. The Wilson function transform: type I. Combining Propositions 14.91 and 14.101 with 
Proposition I4.8( gives the following theorem. 

Theorem 4.12. The Wilson function transform of type I, T : A4 — > Ai, defined by 

(.F/) (A) = / f(x)Mx)dm(x), 
is unitary, and its inverse is given by J-^ 1 = T . 

Proof. First we show that To T is the identity operator on .Mo- For the continuous part of T this 
is Proposition 14.81 therefore we just write down the proof for the discrete part of T. Let g £ .Mo- 
Recall that T> is a finite set, then we obtain from Propositions 14.91 and 14. lUl 

<P\>(x) M 9W4>\(x)ResM(fi) J dm{x) = i ^ #(A)Res M(//) (j)\(x)(t>x'(x)dm(x i 

fo, A'G[0,oo), 
\g(X'), X'eD. 

Combining this with Proposition l4.8l we obtain the desired result. By duality we obtain (F(F/)) (x) = 
f(x) for / £ .Mo- 
Next we show that T is an isometry on .Mo- For simplicity we assume that the measures dm 
and dm are absolutely continuous. From Proposition 14.71 we obtain 

(Ff,Fg)M = lim — / (F/)(x)(F 5 )(x)M(x)dx 

N^oo Z7T Jq 



1 



OO ^00 



lim — — / / /(A) 5 (A')(0 A ,0v>JvM(A)M(A')(iAdA / 



1 

2^ 



o 



/(A / ) 5 (A / )M(A')(iA' = (/,<7)^. 



In case dm and dm have discrete mass points, the proof runs along the same lines. 

So the operator T : .Mo — > -Mo is unitary. Since M.® is dense in M., T extends uniquely to a 
unitary operator on Ad. □ 

Remark 4.13. An interesting special case is the case a = b + c + d— 1. Then (a, 6, c, d) = (a, 6, c, d), 
so the Wilson function transform of type I is then completely self-dual. 

5. The Wilson function transform: type II 

In this section we give another unitary integral transform which has the Wilson function as a 
kernel, by considering the action of the difference operator L on (a dense subspace of) a different 
Hilbert space than in section 0] The method we use is the same method as in section^ therefore 
we omit some details. 

5.1. The Hilbert space H. Let V + C C 5 be the set of parameters a,b,c,d,t satisfying the 
following conditions: 

a = 1 - d, b = c, teR. 
The dual parameters a,b,c,d are still defined by l|2.6|) . and we define the dual parameter t by 

t = l-b- c-t. 

It is easily verified that the assignment (a, b, c, d, t) *— > (a, b, c, d, t) is an involution on V + . Through- 
out this section we assume (a,b,c,d,t) € V + . 
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Let H be the weight given by 

„, . „. , . Y(a ± ix)Y(b ± ix)Y(c ± ix)Y(l — d ± ix) 

H(x) = H(x; a, b, c, d; t) = -i £-i ; \w \ 

sin7r(t ± ix)Y(±2ix) 

Here we use the notation sin(a ± b) = sin(a + b) sin(a — b). The weight is positive for x £ R. 
Observe that a, 6, c, cf; i) = M(cc; a, 6, c, d)/ sin 7r(t±zx), where M is the weight function defined 
in section 0] We assume that H has only simple poles. This imposes conditions on the parameters 
that can be removed afterwards by continuity. Let D + be the infinite discrete set defined by 

V + = {i(t -n) \n £ Z, t - n < 0}. 

We define the measure dh( ) = dh(-;a,b,c,d;t) by 

/s~i poo 
f(x)dh(x) = — f(x)H(x)dx + iC V /(z)ResH(A 
2-7T /n z=a; 

Here C is the normalizing constant 

C = \/sin7r(a + i) sin7r(6 + t) sin7r(c + t) sin-7r(l — d + t). 

Note that C = C. 

For x G T> + we have explicitly, 

i Res H{z) =Y{a + t)Y{a-t)Y{b + t)Y{b-t)Y{c + t)Y{c-t)Y{l-d + t)T(l -d-t) 

z=i(t—k) 

k-t (a - t)k(b - t)k(c - t)k(l - d - t)k 



2t 2 vr 2 (1 - a - t) k (l - b - t) k (l - c - t) k (d - t) k 

and then we see that for (a,b,c,d,t) € V + , the measure dh is positive. For x € supp dh and 
A G supp dh the Wilson function 4>\{x) is real valued. 

We define the Hilbert space TL = H(a, b, c, d; t) to be the Hilbert space consisting of even functions 
that have finite norm with respect to inner product (•, defined by 



(f,g)n = / f(x)g(x)dh(x) 



5.2. The Wronskian. We denote = iC Res H(z), and we define ko to be the smallest integer 

z=i(t— k) 

such that t — ko < 0. For N > and K > ko we define a pairing (•, -)n,k by 

C l' N K 
(f,9)N,K = — / f(x)g(x)M(x)dx + ^ f{i(t - k))g[i(t - k))H k . 

n ^° k=k 

If / and g are real valued functions in Tl, the limits N, K — > oo give the inner product (/, fif}-^. We 
denote limAr_ >00 (/, g)N,K = {f,g)K, assuming that the limit exists. 
For analytic functions / and g we define the Wronskian [/, g] by 

[f,g](k) = [f[i(t - k))g(i(t -k-1))- f{i(t - k - l))g(i(t - - k))H k . 

Proposition 5.1. For even analytic functions f and g we have 

s~t rN+i 

(Lf,g) N)K - (f,Lg) NjK =— / {f(x)g(x - i) - f(x - i)g(x)} B(x)H(x)dx - [f, g](K). 

to Jn 

Proof. We follow the proof of Proposition 14.21 For even functions / and g we have 

(f,g)N,K = -rf f{x)g(x)H{x)dx, 
where Cn k is a contour in the complex plane defined as follows: 
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• Cn,k starts at x = — N and ends at x = N, 

• Cn,k is invariant under reflection in the origin, 

• Cn,k separates the sequence of poles i(a + n), n G Z>o, from the sequence — i{a + m), 
m G Z>o, and similarly for poles of H(x) corresponding to b, c, 1 — d, 

• Cn,k separates the sequence i(t — n), n € Z, n < K , from the sequence i(t — n), n G Z, 
n > if + 1, 

• Cjv,if separates the sequence of poles i(t — n),n< K, from the sequence —i(t — m),m<K, 
Using B{—x + i)H(x — i) = B(x)H(x), we have 

{Lf, 9)n,k ~ {f, Lg)N,K = J- \ [ ~ ( J (f( x ~ i)9(x) ~ f(x)g(x - ij) B(x)H(x)dx. 

A7T \JC N , K Jc N , K +iJ V ' 

Now we make a counterclockwise oriented closed contour by connecting Cn,k and Cm,k + £ at the 
end points, then the integrand I{x) has two poles inside the closed contour; at x = i(t — K) and 
x = -i{t- K -I). So, 

i (1-jU /w * = s (/_r-/D '^t**uv» 

In the same way as in the proof Lemma 14,11 we can show that the first integral on the right hand 
side is equal to the second integral with opposite sign. For the residues we have 

^ Res I(x) = -hf,g}(K), 

Z x=i(t—K) £ 

and, since I{x) = —I(—x + i), 
iC 

— Res I(x) = - lim (x + Ut - K - l))l(-x + i) 

2 x=-i(t-K-l) x-*~i(t-K-l) y ' 

= lim (y-i{t-K))l{y)= Res I{y) = -\[f \ g](K) . 

y->i(t-K) y=i{t-K) 2 

This gives the desired result. □ 
From H{x) = M(x)/sin7r(i ± ix) and ()4.1j) we find, for !/£K, 

H(x + iy) = (x 2a+2b+2c - 2d - 1 e- i7rx ), x -» oo. (5.1) 
This gives the following for the Wilson functions. 
Proposition 5.2. For A ^ X' 

((f>\,(f>\i)K = A , 2 _ A2 • 

Proof. Since the Wilson function is an eigenfunction of L for eigenvalue a 2 + A 2 , the result follows 
from Proposition l5.21 because the integral on the right hand side in Proposition 15 . 21 is equal to zero. 
Indeed, using (|5,1|) and the asymptotic behaviour of (f>\, which follows from Proposition 14.41 and 
(|4.3|) . and applying dominated convergence, we obtain 

i-N+i 

lim / {4>\(x)4>\i(x — i) — <j)\(x — i)4>\>(x)} B(x)H(x)dx = 0. 

A'- *oc J N 

□ 

To find the the asymptotic behaviour of the Wronskian [<f)\, (fty](K) for K — ► oo, we need the 
asymptotic behaviour of Hk, B(i(t — K)) and <fi\(i(t — K)). First we expand in functions with 
nice asymptotic behaviour, as in Proposition 14.41 
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Proposition 5.3. For x = i(t — k) G T> + 

<j> x {x) = d(X)Q x (k) + d(-A)e_ A (fe), 

(-l) fc r(l -5-c-iA-t + fc) / a + iA,6 + iA,c + iA,l-d + *A A 
" A ~~ 7T r(6 + c + iA-t + fc) 4 3 \Ji + c + iA-t + fc,& + c + tA + i- k,l + 2i\ ' J ' 

j . r(— 2iX) sin 7r(t — zA) 

"(A) — 



T(a - iX)T(b - iX)T(c - iX)T(l -d-iX) 
Proof. This follows from Proposition 14.41 and Euler reflection formula; 

I . _^uY(l — b — c — iX + ix) sin7r(6 + c + t + iX) 



T(b + c + iX±ix) tt T(b + c + iX + ix) 

for x = i(t - k), k G Z. □ 
Remark 5.4. Observe that 

- =g(A). 
d(A)d(A) 

For fc — > oo we find from the explicit expression for Q\ and Q4.2JI . for y G Z, 
©a(^ + y) = (zll!^-a-6-c-2a (\ + J_ (d _ a _ 6 _ c _ 2iA )(l -2t + 2y) +o(±j). (5.2) 



7T 

Furthermore, for k — > oo, 

it _ 27f2 i, 2a+26+2c-2d-l A , / 1 

*(*(*-*)) = *(l + (l 

We can find the Wronskian <p X '](K), for K — > oo, from the four Wronskians [0±a, ®±A'](-^0- 
Lemma 5.5. For i"T — > oo, 

[ex,e x ,}(K) = ^(x>-x)K-^^ 

Proof. We have 

[e A , e v ](K) = [@ x (K)@ y (K + 1) - e A (K + i)e v (K)}B(i(i - 

The lemma follows from this expression using the asymptotic behaviour of A (iT), B(i(t — K)) and 

ff*. ' n 

5.3. Continuous spectrum. In this subsection we assume A, A' > 0. 
Proposition 5.6. Let f be an even continuous function, satisfying 

' , X a+b+c-d-±-e e -27r\^ A — > OO, e > 0, 

0(A 5 ), A |0, 5>0. 



/(A) 
Then 



Jim f°° f(X){<t> x ,cf> x ,) N dX- i(A > 
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Proof. From Proposition 15. HI we find, 

6,«6{-l,l} 

and then we obtain from Lemma 15.51 and Proposition 15.21 

i „ K -2i(e\+tX) / / i 

e,€6{-l,l} V 

We multiply both sides with an arbitrary function /(A), and we integrate over A from to oo. The 
function / must satisfy certain conditions that we determine later on. Letting K — > oo then gives 

poo 

lim / f{\){4>\,4>\')Kd\ = 

K^ooJ Q 

lim 1 /°°/(A)(Vi(A)cos(2[A + A / ]ln(K)) + V 2 (A) sin (2[A + A'] ln(ET)) 

K"->oo G J L 



where 



/, , , sin(2[A- A'lln(K)) 1 

+ ^s(A) cos (2[A - A'] ln(K)) + ^(A) ^_ y } 

^(a) = ^ (J(A)d(AO + d(-A)J(-y; 

^( A ) = (d(A)J(-A') - d(-A)d(V; 

^ 4 (A) = -i(d(A)d(-A') + d(-A)d(A') 



riA 



^3 has a removable singularity at A = A'. From the Riemann-Lebesgue lemma we find that the 
terms with fa, i = 1, 2, 3, vanish, provided that G i^ 1 (0, oo). Using 1|4.4|) for the term with ^4, 
we obtain 



00 



lim / /(A)(<Aa,<M^A = 5^ 4 (A')/(A') = ^d(A'M-A')/(A') - ^ 







Using, for i = 1, 2, 3, 4, 

|*(A)| _ A ^00, 

(OtA- 1 ), A|0, 

we see that if / satisfies the conditions given in the proposition, then ftpi € L 1 (0, 00). □ 

Let Ho be the dense subspace of Ti defined by 

Ho = {/ <E H I /(«(t - fc)) = for k » 0} . 

For / € Ho w e define a linear operator Q by 

(0/)(A) = / f(x)Mx)dh(x). 



We call £/ the Wilson function transform of type II. We denote the continuous part of the above 
integral by Q c f, i-e. 

/~1 /'OO 

(Qcf)W = 7T / f(x)4> x (x)H(x)dx. 



2tt„ ., 

From the asymptotic behaviour of 4>\ and H{x) we find that both Qf and £/ c / are well defined if 

feH . 
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Proposition 5.7. Let g be a continuous function satisfying 



5(A) 



O A 5 " 1 ), A |0, 5>0, 



then (F(f c g))(\)=g(X). 

Proof. The proof runs along the same lines as the proof of Proposition 14.81 □ 

Note that if g £ 7Y , then g = o{\ d ~' a - b - t e 2 ' nX ) for A -> oo. So if g £ H , then g satisfies the 
conditions of Proposition 15.71 

5.4. Discrete spectrum. In this subsection we assume A £ P + . In this case d(—X) = and 
3A < 0, and then it follows from Proposition 15.31 and the asymptotic behaviour 1)5. 2 Jl of Q\, that 
4>x £ TL. The following two propositions are proved in a similar way as Propositions 14.91 and 14.101 
therefore we omit the proofs here. 

Proposition 5.8. For A £ V + , \' £ supp(dh), A ^ X' , we have 

(</>A,<M« = o. 

Proposition 5.9. For X £ T> + , 



(<t>\A\)n = (iCR<*H{X')\ 



5.5. The Wilson function transform: type II. In the same way as in the proof of Theorem 
14.121 we can combine Propositions 15.81 and 15.91 with Proposition 15.71 to find that the operator 
Q : Tto Ho is a unitary operator with inverse Q. Since 7^o is dense in 7i, Q extends uniquely to 
a unitary operator on H. 

Theorem 5.10. The Wilson function transform of type II, Q : H — > H, defined by 

(Gf)(X) = ! f(x)<f> x (x)dh{x), 



is unitary, and its inverse is given by Q 



-i 



Remark 5.11. For a = b + c + d— 1, t = ^(1 — b — c — s), s £ R, the Wilson function transform of 
type II is completely self-dual. 



6. Explicit transformations 

In this section we calculate explicitly the Wilson function transforms of certain functions. First 
we give an integral representation for the Wilson function related to Jacobi functions. This leads 
to two explicit transformations in Theorems 16.21 and 16.51 Then, in Theorem 16 .71 we show that the 
Wilson function transform of type I maps an orthogonal basis of polynomials in the Hilbert space 
M to itself, with dual parameters. 

6.1. Transformations related to Jacobi functions. The Jacobi functions, see ^21) are defined 
by 

(eft, ^ P (\{a + (3 + l-iX),\(a + [3 + l + iX) \ 
ip\ ; (x) = 2 i ? il 2 a + ^ ;-xj. (6.1) 
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For x > 1 we use here the unique one valued analytic continuation of the 2-^i-function. The Jacobi 
functions are the kernel in an integral transform pair, called the Jacobi transform, given by 

/>oo 

(Jf)W= I f(z)^' \x)A a , (x)dx 

(6.2) 







where a > -1, f3 € R U iR, 

A Qi/3 (x) = x Q (l + ^, 

and dv{\) is the measure given by 



— / g{X)dv{\) = — \ g(\)\c at p(\)\- 2 d\-i^g(\)Res(c a ,p(v)c at p(-fi)) 
J k Jo XeS t* 

2~ iX r(a + l)r(iA) 



rQ(a + p + 1 + »A))r(£(a - /3 + 1 + »A)) ' 

£ = - a - 1 - 2j) | j e Z> , |/9| - a - 1 - 2j > o}. 

The Jacobi polynomials Pb a '^\x) are orthogonal polynomials on the interval [—1,1] with respect 
to the measure (1 — x) a {\ + x)@dx. They have an explicit expression as 2^1-series, see pQ, [5]. In 
this paper we need Jacobi polynomials of argument (1 — x)/(l + x). Using Pfaff's transformation, 
we find that the explicit expression for these polynomials is 

The orthogonality relation for these polynomials reads, for 3?(a),5?(/3) > —1, 

r ^ (fs) (its) + *^ ib = 



1 T(n + a + l)r(n + ^ + 1) 



(6.4) 



2n + a + (3 + l n\T{n + a + 13 + 1) 

The following proposition gives a representation of a very-well poised 7-Fg-series as an integral 
over a product of two 2-^1-series. 

Proposition 6.1. For (3, fi 6 C, R(a) > 0, £(7 + 5 ± p) > 0, 

\ 2a / \ 2a / 

_ T(2a)r(2a + 5 + p + t)T(J + 7 db p)r(£ - 7 + P) 
T{a + 5 + p±p)T(a + 5 + j±(3) 
x VF(2a + 5 + ^-1 + p;a + j + p,,a + j-p,,a + j3 + p, a - j3 + p, 5 + 7 + p). 

(6.5) 

Proof. We start with the following formula, for !ft(a) > 0, Ui(7 + <5 ± p) > 0, 
^a-ln „A7+<5-a-/3-i_z7V a + / 3 + / 9 ' a + / 3 -/ 5 . V 







V 20 (6.6) 

r(2a)r( 7 + j + p)r( 7 + g - P ) 

"T(a - /3 + 7 + 5)r(a + (3 + 7 + <5) ' 
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To prove this identity we transform the 2-^i-series by Pfaff's transformation ^ Thm.2.2.5], then we 
obtain a 2i ? i-series that converges uniformly on [0, 1] for < 0. We interchange summation and 
integration, then the result follows from using the Beta-integral and Gauss's summation formula p] 
Thm.2.2.2]. The condition on p can be removed using the symmetry in p and —p, and continuity 
in p. 

Next we write the integral in the theorem as 

/= f 1 y 2a-l {l _ y) S^-^2a-l 



J 

where we substituted u » y/(l — y). By [IJ §2.10(3)] the first 2-F\-function can be expanded in 
terms of 2i ? i-functions of argument 1 — y 



a + p + 7, a + p - 7 y 



2a y - 1 

r(2a)r(-2 7 ) M u 1 ^a + 7 + ^,a + 7 -^ 



(l-y) W7 2J Fi ' 1 ";l-y (6.7) 



r(a + /x - 7 )r(a - /i - 7) V 1 + 2 7 

+ r(a + / u + 7)r(a-/, + 7) 1 W V 1 - 2 7 

Observe that the second term is equal to the first term with 7 replace by —7. The integral / splits 
according to this as / = / 7 + /_ 7 . We use formula Q6.6JI to evaluate J 7 . Interchanging summation 
and integration, which is allowed for 3?(a) < 5 since then the 2-^1-series converges uniformly on 
[0,1], and using (|6.6|) leads to 

j. = r(2a) 2 r(-2 7 )r(^ + 7 + P )v(s + 7 - g) 

7 r(a-7 + M )r(a-7-Ai)r(a + /? + 7 + <5)r(a-/3 + 7 + 5) (g ^ 

/a + 7 + /i,a + 7-/i,5 + 7 + /),5 + 7- /) x 
X43 V 27 + l,a + /3 + 7 + 5,a-/3 + 7 + -5 ' 

Then I = I 7 + i_ 7 is the sum of two balanced 4-F3-functions, which by [2J §4.4(4)] can be written 
as the very-well poised yi^-series given in the proposition. 

Note that the expression in (|6.8|) is an analytic function in a for K(a) > 0. The integrand of the 
integral / is analytic in a for 3fJ(a) > 0, and continuous in y. So differentation with respect to a 
and integration with respect to y can be interchanged, see e.g. |22l §4.2]. We see that the integral 
/ is an analytic function in a for 5ft(a) > 0, and therefore the condition 3?(a) < \ can be removed 
by analytic continuation. □ 

Both 2-Fi-series in Proposition 16. II can be considered as Jacobi functions. Using the substitutions 

ai-> -(o + l- d), (3 i— > — (c — b) , 7 i— > -(a + d — 1), 5 i— ► -(6 + c), p^ix, p\-^i\, (6.9) 

and the definition Q3.2JI of the Wilson functions, the right hand side of (|6.5I) can be written as 

r(l + a - d) 2 T{a ± iA)r(l -d±i\) <p x (x; a, b, c, d). 

So Proposition 16.11 gives a representation of the Wilson function as the Jacobi function transform 
of a Jacobi function. From the inverse Jacobi transform we obtain the following. 
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Theorem 6.2. Let tp(x) = </4t +C 1,c b+1 \u), u>0, and ip(x) = sin-7r(t ± ix)tp(x), then 

(00(A) =C(l + u)i- &+ ^ 2 F 1 (~ a + iX ^~ d + iX -_ u \ (a,b t c,d,t) 6 V + . 

\ 1 + a — a I 

Observe that for u = the first statement gives (J-1)(X) = 1. 

Proof. In Proposition 16 . II we replace p by iX G iR. Then using the definition of the Jacobi functions 
(|6.1|) we can write Proposition 16, II as 

y (l+uT P 2-Fil 2a ;~ U J^2A («) A 2a-l,2/3 («)» = F(\), 

where F(X) denotes the right hand side of 1)6. 5 Jl . Note that from 1)6.7)) we find for u — > oo 

(1 + u)^ 2Fl ^ + M + 7, « + M - 7 . \ „ ^-o-ZJ-*-, + C 2 «— 

where Ci and C2 are independent of u. So for ±7) > we see that this function is an element 
of L 2 ([0, 00), A2 a -i 2p(,u)du). Taking the inverse Jacobi transform, assuming for simplicity that the 
discrete set £ is empty, gives 

(i + «)^-^if° + /4 + 7 ' a + /1 ~ 7 ;- 

\ 2a 



1 

2^ 



F (x)^r^\u) 



T(a + (3 + iX)T(a - (3 + iX) 



T(2a)T(2iX) 



2 



dX. 



We write out F{X) and ip^x l ' 2 ^\t) as hypergeometric functions and use the substitutions given 
in ()6.9|) . without the last substitution, then we obtain 

-c+ix ^ { a + ix,l — d + ix \ 1 f°° ^ fc + iX,c — iX 



This is the first statement in the theorem with dual parameters, and x and A interchanged. Writing 
the above integral as a contour integral (with the contour as in the proof of Proposition 14.21 with 
N — > 00), and using analytic continuation in the parameters a, b, c, d, the identity can be extended 
to the case (a, b, c, d) € V. Deforming the contour again to the real line might add a finite number 
of discrete mass points. 

The second statement follows from the first. We use H{x) = M(x)/sin7r(t±ix) and we use that 
the function sin7r(i ± ix)ip(x) vanishes on the infinite discrete set T> + . □ 

We substitute \i = —a — 7 — n, n £ Z>o, in Proposition 16.11 then the first 2-^i-series on the left 
hand side of (|6.5[) terminates and can be written as a Jacobi polynomial by 1)6.3)) . Writing the right 
hand side as J 7 + J_ 7 , with I 7 as in (|6.8j) . we see that J_ 7 = because of the factor r(a + 7 + /x) _1 . 
Now we obtain, for K(a) > and 5R(7 + 5 ± p) > 0, 



J°° u 2-i (1 + u) /3-^ 7p (2 Q -i,2 7) n_z_^ ^ + /3 + P, * + " P . _^ 



it du 



_ (-1)"(2 7 + l)nr(2a)r( 7 + g + g)T(7 + 5 - p) 

n\r(a + P + j + 5)T(a-(3 + ~f + 5) 1 ' ' 

/-n, 2a + 27 + n, J + 7 + p, 5 + 7 - p 
X4 3 V 2 7 + l,a + /3 + 7 + ( J ! a- y 5 + 7 + <5 ' 
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This is Koornwinder's formula |131 (3.3)] stating that Jacobi polynomials are mapped onto Wilson 
polynomials by the Jacobi function transform. We use this formula to expand the integral in 
Proposition 16. II in terms of terminating ^^i^-functions. This leads to an expansion formula for the 
Wilson function. 



Proposition 6.3. Let f,g G C satisfy > 0, »(/ + g) > 0, 3ft(a + g) > 0, »(1 - d + 5) > 0. 

T/ien 

, / , ,s / u „ f-n,a + f + g-d + n,f + i\. 

<px{x;a,b,c,d) = > C n {x,\)4F 3 ~ _ 

V / + 9,b + f,c+ f 

-n, a — d + f + g + n, g + a, g + 1 — d 




f + g, g + b + c + ix, g + b + c — ix ' 
where 

T(a-d+f + g + l)T(g + a)T(g + 1 - d)T(f ± iX) 



C n (x, A) 



r(l + a - d)T(f + g)T(f + b)T(f + c)T(g + b + c± ix)T(a ± i\)T(l -d±i\) 
a + f + g-d + 2n(a + f + g-d) n (f + g) n 



a + f + g - d n! (a + 1 - d) n 

Proof. Let I be the integral in Proposition 16.11 We expand the first 2-^1-function in I in terms of 
Jacobi polynomials of argument (1 — u)/(l + u): 

(1 + uT +^-^ 2Fl (<* + * + 7. « + ^ " 7 . _ \ = £ Cn p^-xM 

\ a* J n=Q \ +uj 

Note that from (|6.7j) we obtain that the left hand behaves for large t as Ciu ri ~ a ~' y+1 + C2?/ ? ~ cr+7+1 , 
where Ci and C% are independent of So for 5R(r/ + cr ± 7) > the function on the left hand side 
is an element of L 2 ((0, 00), u 2 "" 1 (1 + u)~ 2a ~ 2ll ~ 1 du) . We see that the expansion on the right hand 
side converges uniformly for u in compact intervals. The coefficients c n are found using ()6.1Uj) and 
the orthogonality relation for the Jacobi polynomials (|6.4|) . and this gives for 5ft(a) > 0, > — ^, 
and 3?(<r + 77 ± 7) > 0, 

_ r(2a + 2?? + l)r(? ? + (T±7) (2a + 2?? + 2n) (-l) n (2a + 2r?) n 

° n ~ r(2r/ + l)r(a + 7/ + a ± /x) (2a + 2t/) (2a)„, 

/— n, 2a + 2?? + n,a + 77 + 7, a + 77 - 7 
X 4-T3 I ; 1 

\ 27] + l,a + fj, + rj + a,a — n + rj + a 

From this expansion we find 

1 = J> jf + uf-«- s -^P^-^ (1^) 2 F X ( a + " + ^ + ^ - P ; -„) 

Using (|6.1(JI) again gives the following expansion for I, for »(a) > 0, K(r/) > -|, 3f?(cr + 7/ ± 7) > 0, 
5R(<5 - a + r/ + 1 ± p) > 0, 

_ r(2a)r(2a + 2r? + l)r(r/ + a ± 7 )r(r/ + a + 
r(2r7 + l)r(a + 77 + a ± p)T(a + ?7 + <5-<t + 1±/?) 

2a + 2r] + 2n (2a + 2r]) n (2r] + l) n /-n, 2a + 2t/ + n, a + 77 + 7, cr + 77 - 7 



00 

x 

n=0 



' 2a + 2r/ n! (2a) n 4 3 \27/ + l,a + /i + ?7 + (T, a — /i + 77 + 0" 

/— n, 2a + 277 + n, 5 - a + 1 + 77 + p, 5 - a + 1 + 77 - p ^ 
X 4 3 V 277 + l,a + /3 + ?7 + 5-a + l,a-/3 + ?7 + (5-(j + l ' 
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The proposition then follows from Proposition l6.ll the substitutions g — a + 1 + 5 \— > /, g + a — 5 i— > 5 
and □ 

Corollary 6.4. For &(/) > 

<f)\(x;a,b,c,d) = 

/ u I? f-nJ + iXJ-iX \ ( -n,f-a,f + d-l \ 
^ V b + f,c + f J yf + c-c-ix,f + c-c + ix J 

where 

(/ + c- c + ix) n {f + c- c-ix) n T(f + iX) 

C n (X,A) = = = . 

n\ r(l + a - d + n)T(f + b)T(f + c)T(a + iX)F(l — d±iX) 

Proof. We use Whipple's transformation ^ Thm.3.3.3] to write the second 4i ? 3-function in Propo- 
sition 16,31 as 

(/ + c-c + ix) n (f + c-c- ix) n ( -n, a - d + f + g + n, f - a, f + d - 1 \ 
(g + b + c + ix) n (g + b + c + ix) n 4 3 ^ f + g, f + c - c - ix, f + c - c + ix ' J' 

Then the corollary follows from letting g — > oo, using (|4.2|) for the T-functions. □ 

We consider one of the 4-F3-functions in Proposition 16.31 as a Wilson polynomial, then we obtain 
the following Wilson function transform of a Wilson polynomial. 

Theorem 6.5. Let f,g€M. such that g+a, g+b, g+c, g+l — d has positive real part. Furthermore, 
let (p n and ip n be the function defined by 

ipn{x) = T(g ±ix)R n (x; f,b,c,g), ip n (x) = sin7r(t ± ix)<p n (x). 

Then 

CF, n \(\\ r-i(r,h \(\\ r(ff + o)r(g + b)T(g + c)T(g + 1 - d) (g + b) n (g + c) n 

(Fip n )(X) = C {Gipn){X) = — = — — / j , t \ (f . \ 

T{g + b + c±i\) U+b)n{f+c) n 

f-n,b + c + f + g + n-l,g + l-d,g + a 
x 4F3 ~ ~ ; 1 

V f + g,g + b + c + iX,g + b + c-i\ 

Here we assume for T that (a, b, c, d) £ V, and for Q that (a, b, c, d, t) G V + . 

Proof. We recognize the first 4-F3-function in Proposition 16.31 as the Wilson polynomial R n (X) = 
R n (X; f, b, c, g). We multiply by R n (X) and integrate against the orthogonality measure d[i(X; f, 6, c, g). 
Using the orthogonality relation (|2.3|) then gives 

-J o R n (X)M^a,b,c,d)^ U L__U = 

r(g + q)r( g + b)r(g + c)T(g + 1 - d) (g + b) n {g + c) n 
T(g + b + c + ix) (f + b) n (f + c) n 

( -n,a- d + f + g + n,g + 1 - d,g + a \ 
V f + g,g + b + c + ix,g + b + c — ix I 

The theorem follows from this by replacing the parameters a,b,c,d by their dual parameters, and 
interchanging A and x. 

Note that the integral can again be written as a contour integral. Then by analytic continuation 
the result remains true when g € R and the pairwise sum of a, b, c, d, g has positive real part. □ 

From Corollary 16.41 a similar theorem can be derived for the Wilson function transform of a 
continuous dual Hahn polynomial, which is a limit case of the Wilson polynomial, see pQ, [H]. We 
leave this to the reader. 
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6.2. The Wilson function transform of type I of an orthogonal system. Koornwinder's 
formula (|6.10j) shows that the Jacobi function transform maps a (bi) orthogonal system of polyno- 
mials (the Jacobi polynomials) in L 2 ((0, oo), A ai p(t)dt) , onto an orthogonal system of polynomials 
(the Wilson polynomials) in L 2 ((0, oo), du) . We show that there exists a similar formula for the 
Wilson function transform of type I. 

The Wilson polynomials R n {x) = R n (x;a,b,c,l — d) form an orthogonal basis for the Hilbert 
space Ai. We show that the Wilson function transform of type I maps a Wilson polynomial R n 
onto itself, with dual parameters. The proof is based on the following proposition. 

Proposition 6.6. Let L be the difference operator as in Proposition \H. 1\ then 

(LR n )(x) = C n R n+1 (x) + (C n + D n )R n (x) + D n R n ^(x), 

where 

(n + a + b + c — d){n + a + b)(n + a + c)(n + a + 1 — d) 



C n 



(2n + a + b + c - d)(2n + a + b + c - d + I) 
n(n + b + c — l)(n + b — d) (n + c — d) 



(2n + a + b + c - d)(2n + a + b + c - d - 1) 

Proof. Applying the difference operator L on R n gives 

{LR n )(x) = B(-x)R n (x + [A(-x) + A(x)] R n (x) + B(x)R n (x - i) 

= [n(n + a + b + c-d) + A(x) + A(-x) + B{x) + B(-x)] R n (x). 

Here we use the difference equation (j2.4j) for the Wilson polynomials R n (x; a, b, c, 1 — d). From the 
explicit expressions for A(x) and B(x) we obtain 

A(x) + B(x) = (a + **)(fe + *'0(c + *'0 > 

and this gives 

(a + ix)(b + ix)(c + ix) — (a — ix)(b — ix)(c — ix) 



2ix 

ab + ac + bc — x 2 . 



A(x) + A{-x) + B(x) + B(-x) 



Then, using the three-term recurrence relation (|2.1|) for the Wilson polynomials, we find 

(LR n )(x) = [n(n + a + b + c — d) + ab + ac + be + a 2 ] R n (x) — (a 2 + x 2 )R n (x) 
=C n R n+ i(x) + D n R n -i(x) 

+ [n(n + a + b + c — d) + ab + ac + be + a 2 — C n — D n ]R n (x). 

Note that we use here C n and D n as in ()2.1|) . but with d replaced by 1 — d. A long, but straight- 
forward calculation shows that 

C n + D n + C n + D n = n(n + a + b + c — d) + a 2 + ab + ac + be, 

and from this the proposition follows. □ 

Theorem 6.7. For (a,b,c,d) G V, the Wilson function transform of type I of the Wilson polynomial 
Rn(x) = R n (x; a, b, c, 1 — d) is given by 

= ( - ir (r^)^ (A) - 
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Proof. We use (J r R n )(\) = limjv-«x)(0A> Rn)N-> where (-,-}n is the pairing defined in sectional 
From Proposition 14.21 we find, 

lim {Lcj)\,R n ) N - {(j)\,LR n ) N = lim [<px, Rn](N). 

N-+00 N^oo 

Using (|4.1|) . (|4.3j) . R n (x) = 0{x 2n ), and applying dominated convergence, we obtain for the Wron- 
skian 

lim [<j> x ,R n ](N) = U. 

The Wilson function (p x is an eigenfunction of L for eigenvalue a 2 + A 2 , so we have 

{d 2 + \ 2 ){FR n ){\) = lim (L(j> x ,R n ) N = lim (<f) X ,LR n ) N = (F(LR n ))(\). 



Then from Proposition 16.61 and linearity of J-, it follows that the function TR n satisfies the recur- 
rence relation 

(a 2 + A 2 ) y n (X) = C n y n+ i{\) + (C n + D n ) y n (X) + D n y n _i(A). (6.11) 

From R^(x) = 0, we obtain (^E_i)(A) = 0. By Theorem l6~2l we have (TR )(\) = (J r l)(A) = 1, 
so we find from the recurrence relation that (TR n ){\) = P n {\). Here -P n (A) is a polynomial in 
A 2 of degree n satisfying the three-term recurrence relation (|6.11j) . with initial values P_i(A) = 
and -Po(^) = 1- From the recurrence relation for R n (X) = R n (\;a,b, 5,1 — d) it follows that 

P nW = ("If (T^^n(A). □ 

Remark 6.8. Theorem 16.71 gives in fact a new proof for the unitarity of JF, since an orthogonal basis 
in M. is mapped to an orthogonal basis in Ai with the same norm. Note that Theorem 14. 121 is not 
used to proof Theorem 16.71 

Theorem 16.71 can be considered as the q = 1 analogue of 21, Prop.4.1], which is used by Stok- 
man to obtain an expansion of the Askey- Wilson function in Askey- Wilson polynomials, see |21[ 
Thm.4.2]. So one might expect that from Theorem 16. 71 a similar expansion can be proved for the 
Wilson functions. However, if we formally expand 

oo 

4>x(x; a, b, c, d) = d n (x, A) R n (x; a, b, c, 1 — d)R n (X; a, b, c, 1 — d), 

n=0 

and we calculate d n (x,\) using Theorem 16.71 and the orthogonality relation for the Wilson poly- 
nomials, we find an expansion that in general does not converge, but it does converge if x G T>, 
or A G T>. The convergence of the expansion formula for the Askey- Wilson functions is due to a 
Gaussian factor g m ( m + 1 )/ 2 ) < q < 1. In the limit q — > 1 this factor disappears. 

From Theorem 16 . 71 we can also find the Wilson function transform of type II of a Wilson poly- 
nomial. 

Theorem 6.9. For (a,b,c,d,t) G V + , the Wilson function transform of type II of the Wilson 
polynomial R n (x) = R n (x; a, b, c, 1 — d) is given by 

(GR n )(X) = {-IT, {b + C \ - sinvr^A) A 

(1 + a - d) n ^sin vr(a + t) sin ir(b + t) sin vr(c + t) sin tt(1 -d + t) 

Proof. From Theorem 16.71 we obtain 

(9 sm7r(t±ix)R n )(X) = (-If n C(b + C j," RnW- 

(l + a-d) n 

Taking the inverse of this gives 

(gR n )(x) = (-1)" ~ f" sinvr(t ± ix)R n {x). 

O [0 + C) n 
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Replacing all parameters by their dual and interchanging x and A then gives the statement in the 
theorem. □ 
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